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An  Integrated,  Optimization-Based  Approach 
to  the  Design  and  Control  of  Large  Space  Structures. 


Elijah  Polak.  Karl  S.  Pister,  Robert  L.  Taylor 

PROGRESS  REPORT  for  tha  I^riod  October  1.  1983  to  May  1,  1984. 

The  following  activities  took  place  in  the  report  period. 

1.  In  January  1984,  we  offered  a  three  day  workshop  on  the  fundamentals  of 
semi-infinite  optimization  in  engineering  design.  The  workshop  was  attended  by 
20  participants  from  Air  Force  Navy  and  NASA  laboratories,  some  of  their  con¬ 
tractors  and  a  few  civilians  from  industry.  A  memorandum  [l]  Iwas  prepared  for 
distribution  to  the  participants.  This  memorandum  presents  a  new  theory  deal¬ 
ing  with  the  construction  of  semi-infinite  optimization  algorithms  for  engineer¬ 
ing  design.  A  sequel  to  this  memorandum  is  is  currently  being  prepared.  It  is 
expected  that  the  two  memoranda  will  serve  as  a  basis  for  an  invited  review 
paper  to  be  published  by  SIAM  reviews. 

2.  We  have  examined  various  alternatives  in  the  formulation  of  the  optimal 
integrated  design  of  a  large  space  structure  and  its  control  system.  We  have  set¬ 
tled  on  a  model  problem  which  is  to  serve  as  a  basis  for  our  research.  This  model 
problem  contains  a  number  of  features  that  have  previously  not  been  dealt  with 
in  the  context  of  optimization-based  design  and  it  has  led  to  a  number  of  very 
challenging  and  exciting  new  problems.  Our  model  problem  is  presented  in 
Appendix  1. 

3.  Three  papers  have  been  completed  which  deal  with  various  aspects  of 
optimization-based  control  sjrstem  design. 


In  [2]  we  describe  our  recent  work 
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opUmization  based  control  system  design  zis  well  as  the  progress  that  has  been 
made  on  our  control  system  design  software  package,  DELIGHT.  MIMO. 

In  [3]  we  present  a  survey  of  the  basic  aspects  of  our  methodology  of  design  of 
linear  multi-variable  control  systems  via  semi-infinite  optimization.  Specific 
topics  treated  are  (i)  data-base  and  simulation  requirements,  (ii)  techniques  for 
the  transcription  of  design  specifications  into  semi-infinite  inequalities,  and  (iii) 
lemi-infinite  optimization  algorithms  for  control  system  design. 

In  [4]  We  show  that  wor«!t  case  design  of  control  systems  vrith  both  parametric 
and  unstructured  uncertcxinty  leads  to  a  new  class  of  semi-  infinite  optimization 
problems.  We  present  a  three  phase  seiiii-infinite  optimization  algorithm  based 
on  the  theory  developed  in  [l]  for  solving  this  class  of  problems. 
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APIVNDIX 1:  MODEL  PROBLEM  OP  AN  INTEGRATED  LSS  OPTIMAL  DESIGN. 

(i)  IntroducUoD. 

Aa  we  have  already  stated  in  the  introduction,  we  propose  to  consider  the 
design  of  large  space  structures  which  are  required  to  perform  large  amplitude 
maneuvers  at  the  end  of  which  they  are  required  to  remain  "locked"  on  a  target. 
We  propose  to  deal  with  the  the  pointing 

of  the  LSS  in  two  stages.  In  the  first  stage,  the  control  task  is  to  rapidly 

redirect  the  pointing  direction  of  the  LSS  reference  axis,  e.g,.  the  line-of- 
fight  of  a  telescope  or  antenna,  to  a  desired  target  direction.  At  the  end  of 
the  first  stage,  which  is  to  be  carried  out  by  open  loop  optimal  control,  the 
large  motions  of  the  LSS  induced  by  the  maneuver  must  "quiet  down"  so 
that  the  control  of  the  LSS  can  be  transferred  to  a  linear,  closed  loop  con¬ 
trol  system.  The  task  of  the  latter  is  to  damp  out  the  induced  structural 
vibrations,  to  suppress  disturbance  efiecls  (e.g.,  caused  by  running  cooling 
water  through  pipes),  and,  finally,  to  lock  the  pointing  direction  on  the  tar¬ 
get. 

For  the  purpose  of  obtaining  a  tractable  model  problem  for  our 
research,  we  shall  initially  assume  that  the  LSS  is  a  beam.  We  shall  use  a 
nonlinear  beam  model  for  the  large  motions,  and  we  shall  use  a  linear  beam 
model  to  describe  the  small  displacements  as  a  perturbation  around  the 
equilibrium  rigid  body  configuration. 

In  fully  developing  our  model  problem,  below,  we  shall  first  develop  the 
equations  of  motion  for  a  beam  under  large  and  small  displacement  condi¬ 
tions.  We  shall  then  use  the  resulting  equations  in  transcribing  sample 
design  specifications  into  infinite  systems  of  inequalities. 


(ii)  General  CooaideFations  in  Modeling. 

The  simulation  of  a  large  space  structure,  under  conditions  outlined 
above,  involves  the  analysis  of  a  system  which  is  subjected  to  large  motions 
about  the  center  of  mass  and  to 

large  deformations.  In  order  to  model  the  LSS  response  during  manoeuvers 
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and  to  couple  effectively  with  appropriate  control  and  design  algorithms, 
the  structural  model  must  be  chosen  carefully.  There  are  a  number  of 
aspects  to  this  problem  which  deserve  special  consideration  and  research. 
For  space  structures  which  have  large  motions  with  respect  to  the  center  of 
mass,  the  "rigid  body"  equations  for  angular  momentum  will  contain  time 
varying  "inertia"  parameters.  Accordingly,  an  analysis  procedure  which 
simulates  the  response  in  terms  of  the  "rigid  body"  response  with  super¬ 
posed  large  motions  may  not  be  the  most  efficient  approach.  Furthermore, 
even  in  the  simplified  situation  when  a  small  motion  is  superposed  on  large 
rigid  body  motions,  significant  terms  may  be  omitted  if  the  theory  is  not 
property  deduced  from  a  linearization  of  the  more  general  non-linear 
forms. 

For  the  purpose  of  getting  preliminary  research  results,  we  will  use  a 
large  displacement  three  dimensional  beam  model  which  includes  the 
salient  features  for  modeling  the  response  of  LSS.  A  theory  was  developed 
by  Kirchhoff  and  described,  for  the  static  case,  in  the  treatise  on  elasticity 
by  Love  [Lov.l].  Other  significant  contributions  to  the  theory  of  beams  are 
contained  in  the  references  [Ant.l,  Gre.l,  Nag. 2,  and  Sim.l].  In  the  discus¬ 
sion  which  follows  we  include  the  appropriate  inertia  terms  in  the  theory. 

(lii)  Equations  of  Motion  of  a  Beam;  l^rg e  Displacements. 

A  typical  LSS  manoeuver  is  shown  in  Figure  1  where  a  structure  which 
is  initially  assumed  to  be  in  a  "quiet"  state  is  subjected  to  control  forces 
which  are  to  produce  a  final  state  which  is  also  "quiet". 

During  the  intermediate  times  of  the  manoeuver,  the  motions  in  the  struc¬ 
ture  may  produce  large  relative  displacements,  as  shown  in  Figure  1  for 
time  t'. 

Far  the  simplified  beam  model  of  a  space  structure  this  motion  may  be 
described  by  the  relation  (e.g..  see  Rgure  2) 
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y(t.*)  =  fo(«.*i)  +  f)  ^Ti(t.zi)  (2.1) 

t  B  e 

where  j  describes  the  deformed  position  of  the  beam,  in  terms  of  the  vec¬ 
tor  x=  (>1.22.23)..  where  Z|  c  [0. f]  is  a  parameter  defining  position  along 
the  axis  of  the  beam,  >2,23  are  coordinate  positions  in  the  cross  section  A, 
(i.e..  the  vector  (22.23)  CAclR^),  along  the  directions  defined  by  the  unit 
vectors  T2(f.2i)  and  T3(f.z,),  t  is  time,  and  {Po(f.Zi)  is  the  deformed  posi¬ 
tion  of  the  beam  axis.  We  assume,  for  simplicity,  that  the  tangent  to  the 
deformed  axis  of  the  beam  satisfies  the  first  Frenet  formula 


Ti(f.2,)  =  g7p»o(^Zl) 


(2.2) 


where  - - is  the  derivative  along  the  beam  axis,  and  the  unit  vectors  t.  in 

dzi  * 

the  deformed  beam  are  mutually  orthogonal  (i.e.,  plane  cross  sections 

remain  normal  to  the  beam  axis  during  all  deformations).  Accordingly,  we 

obtain  the  kinematic  theory  of  beams  due  to  KirchhofT  which  is  described  in 

Love's  classical  treatise  [Lov.  l]  and  summarized  below  (for  an  expanded 

description  including  effects  of  shear  deformation  see  [Sim.2]).  The  motion 

for  this  theory  is  completely  characterized  by  fpQ  and  Ti  ,i  =  1,2,3,  as  shown 

in  Figure  3. 

The  balance  of  linear  momentum  for  this  theory  is  given  by  the  relation 


g|Y^(t.*i)  +  q(t,2|)  +  Ui(f,2,)  =p.4(2,)^w(f,z,)  (2.3) 

and  the  angular  momentum  equation  by 

+  57pi>o(«.*i)xf{<.*i)  +  ue(f,2i)  =  /)I(*i)  ^(f.2,)  .  (2.4) 
where  f(f,2|)  and  m(f,2i)  are  the  stress  resultants  and  stress  couples  (as 
shown  in  Figure  4).  <|(f.S|)  is  a  specified  loading  (e.g.,  gravity),  U|(f  ,Zi)  and 
ll2(f,S|)  are  the  control  forces  (e.g..  see  Figure  5  for  types  of  Ui  and  Us 
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quantlties),  p  is  the  mass  density  of  the  beam  material,  is  the  cross 

sectional  area,  I(xi}  is  the  matrix  of  cross  section  moment  of  inertia, 
and  o{t,Zi)  are  the  translational  velocity  and  the  angular  velocity  of 
the  local  triad  Tj  cross  section,  respectively,  and  ^-denotes  material  time 

diflerentiation.  In  the  above,  the  accelerations  are  defined  in  terms  of  the 
rate  of  change  of  velocity  (which  are  the  rates  of  change  of  the  motion) 

through 

»(<.*»)  =  (2.5) 


and 


=  UXT,(f,2,)  (2.6) 

It  should  be  noted  that  this  relationship  may  be  parameterized  through  the 
use  of  Euler  angles  or,  to  avoid  numerical  singularities,  by  using  quatern¬ 
ions  (see  discussion  on  rigid  body  motion  in  section  on  linear  beam  models). 

As  noted  previously  the  kinematic  response  of  the  beam  is  described 
by  f  0  •u'd  the  t*  unit  vectors.  The  t*  may  be  defined  in  terms  of  the  direc¬ 
tion  cosines  with  respect  to  an  inertial  coordinate  frame  described  by  the 
fixed  unit  vectors  e^,  i  =  1,2,3  of  the  inertial  frame  of  reference.  Accord¬ 
ingly. 


Cl 

II 

N 

ez 

es 

(2.7) 


where  A(f  ,Z])  is  the  matrix  of  direction  cosines.  The  rate  of  twist  <i  and 
the  changes  in  curvature  Cg  and  xg  of  the  beam,  which  describe  the  defor¬ 
mations,  may  be  deduced  by  differentiating  the  Ti's  along  the  beam  axis. 

Thus, 


A(<.X|) 

OZi 


A’'(t.*j)T<(f.s,) 


(2.B) 
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where  the  the  twist  and  changes  in  curvature  are  obtained  from 

*(^•*1)  =  =  -«3  0  Ki  (2.9) 

[/cz  -/c,  0 

The  similarity  between  rigid  body  mechanics  emd  the  beam  theory  con¬ 
sidered  here  is  very  evident.  Indeed,  the  beam  axis  parameter  Zj  and  time 
t  lead  to  a  direct  analogy  between  quantities  relating  rales  of  change  along 
*i  and  time  rates  of  change  for  the  unit  vectors  Tj.  Accordingly  methods 
which  accurately  solve  for  the  angular  velocity  in  the  rigid  t  y  motion 
problem  may  also  be  useful  in  solving  for  the  beam  rotations. 

For  the  kinematics  defined  by  (2.1).  the  strains  7(1.  z,)  ar  d  led  by 
the  relation 

7(f.2,)  ^  A^(f.z,)  g|pPo(f,2,)  -  T,(<.z,)  (2.10) 

where,  due  to  (2.2).  the  only  nonzero  component  will  be  the  axial  strain. 

Appropriate  constitutive  equations  to  interrelate  the  beam  stress 
resultants  and  stress  couples  with  the  strains  and  deformations  may  be 

expressed  by 


f  =  f(7.x.z,.f) 

(2.11a) 

m  =  m(7,c.zi,f) 

(2.11b) 

where  f  and  m  are  functions  of  the  arguments  for  elastic  materials  or, 
alternatively,  they  are  functionals  of  the  arguments  for  inelastic  (e.g., 
viscoelastic)  materials  which  can  characterize  dissipation  of  energy  (i.e., 
damping).  As  a  special  case  one  may  assume  linear  elastic  relations 
between  the  stress  resultant  and  beam  axis  strain 


f  =  By. 


(2.12a) 
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f 


and  similarly  for  the  stress-couples  and  changes  in  curvatures  and  twist 

m  =  Dk  (2.12b) 

where  B  and  D  are  matrices  for  the  beam  which  are  constant  in  time.. 

The  response  of  the  beam  may  be  determined  by  solving  the 
difTerential  equations  defined  above  together  with  appropriate  boundary 
and  initial  conditions.  The  initial  conditions  for  a  typical  manoeuver  may  be 
specified  by  defining  the  state  vector  x(t,Zi)  at  an  initial  time  f°.  The  state 
vector  is  given  by  the  position  quantities  ipoit  ,Zi)  and  Tj(<  ,  z, ) ,  i  =  1.2.3 
and  the  rates  of  change  of  the  position  t  ,  z  j )/  and  of  the  unit  vec¬ 
tors  flT,( f  ,  z , )/  flf :  that  is 


*  A  T  -r  T 

-  (^O-Ti.  2,  3. 


9T3  , 
dt  ■ 


(2.13) 


The  boundary  conditions  for  a  typical  LSS  simulation  will  consist  of  equa¬ 
tions  which  specify  the  stress  resultants  and  stress  couples  at  the  ends  of 
the  beam.  The  boundary  conditions  may  define,  ‘or  example,  beams  with 
attached  rigid  end  masses  (we  refer  to  this  model  as  a  "dumb  bell"  struc¬ 
ture).  'n^e  appropriate  boundary  conditions  for  this  case  specify  end  loads 
for  the  beam  from  the  inertial  behavior  of  each  rigid  end  mass. 

In  general,  the  response  of  the  beam,  at  any  time,  may  be  specified  in 
terms  of  the  kinematical  quantities  substituting  the  kinemati- 

cal  measures  into  the  constitutive  equations,  and  then  substituting  the 
resulting  constitutive  equations  into  the  momentum  balance  equations.  If 
we  select  appropriate,  finite  dimensional  design  parameters,  elR"*,  for 
the  beam  (e.g.,  cross  sectional  area,  A,  and  moment  of  inertias,  I,  etc.)  and 
PkCR"*  for  the  actuator  locations,  we  obtain  a  first  order  difTerential  equa¬ 
tion  of  the  form 


*(f.s,)  =  h(i<f  .s,).u,(f.s,).uz(f.z,),pb  .p*,)  . 


(2.14) 
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wbere  h  is  a  dilTerential  operator  in  Zj.  The  control  quantities  Uj( <  ,  Zi )  and 
U8(t  1*1)  are  to  be  used  in  positioning  the  space  structure  into  a  particular 
attitude. 


(iv)  Equations  of  Motion  of  a  Beam;  Sknall  Displacements. 

In  the  second  phase  of  the  manoeuver  the  structural  response  may  be 
modele(^  by  linearizing  the  motion  of  the  LSS  about  a  rigid  body  response. 
Accordingly,  the  beam  model  cjin  be  deduced  by  describing  the  small  dis¬ 
placements  as  a  linearization  of  the  large  displacements  relative  to  a  single 
point,  the  center  of  mass,  in  the  beam  which  describes  the  rigid  body  equili¬ 
brium  conflguration.  In  the  following  description  we  let  a,(f),  i  =  1.2.3, 
denote  the  unit  vectors  of  the  moving  frame  attached  to  the  center  of  mass 
of  the  beam  structure  (i.e.,  Oiit)  =  Ti(t,Zi„/)  i  e..  z,  „/  is  the  position  for 
the  center  of  mass)  and.  as  before,  let  Cj,  i  =  1,2,3.  be  the  unit  vectors  of 
the  inertial  frame  of  reference.  In  addition  the  beam  axis  tangent  vector, 
ai(f).  defines  the  reference  axis  which  is  required  to  point  at  the  target. 

The  solution  of  this  phase  of  the  problem  requires  two  parts:  (a)  the 
computation  of  the  rigid  body  response  (which  may  still  involve  large  posi¬ 
tion  and  angle  changes  if  the  tracking  times  are  of  long  duration):  and  (b) 
the  computation  of  the  beam  response  about  the  rigid  body  equilibrium 
state. 


Tlie  rigid  body  motion  of  a  beam  about  its  center  of  mass  is  described 
with  respect  to  the  moving  coordinate  system  ai{f ),  by  two  equations.  The 
translational  equation  of  motion  is 


_  <11(0  +  «.(ui(f).Pi,) 

^0 - jf - 


-  6>(f)xv(t)  , 


(2.15a) 


and  the  rotational  equation  of  motion  : 
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6(0  =  J  Mqzi-l  +  it(ue(0.Pb)  -  «X(J«) 


(2.15b) 


where  U  is  the  total  mass  of  the  beam,  J  denotes  the  inertia  matrix,  qj  ,  Uj 
and  4; ,  Ug  are  applied  and  control  forces  and  couples,  respectively. 
■i(*ii(0tPh)  “id  “aCuiCO.UaCO.Pi)  denote  the  resultant  control  force  and 
couple  on  the  rigid  body,  respectively,  u  and  v  represent  the  rigid  body 
angular  and  linear  velocities,  respectively. 

A  rate  equation  describing  the  evolution  of  the  moving  frame  with 
respect  to  the  Axed  inertial  frame  may  be  written  in  the  form 

X  =  |-nx.  (2.16) 

where  x  denotes  a  quaternion  and  0  is  a  skew-symmetric  matrix  whose 
coelTicients  are  composed  of  the  components  of  the  angular  velocity,  w*  as 
follows 


0 

cJa 

-  Wg 

-  CJ3 

0 

Wg 

Ug 

-  w, 

0 

Ws 

*  Wg 

‘U3 

0 

The  use  of  a  quaternion  avoids  singularities  which  are  often  encountered  in 
alternative  methods  (e.g..  use  of  Euler  angles). 

As  noted  above,  in  the  second  part  of  the  analysis  we  consider  a  lineari¬ 
zation  of  the  large  displacement  beam  theory  about  the  rigid  body  response 
just  described.  In  this  analysis  we  let  iJj^(f,C|),  i  =  1,2,3,  denote  the 
linearized  beam  displacements  (with  respect  to  the  rigid  body  motions) 
along  the  principal  axes  The  beam  is  assumed  to  undergo  both  axial 

and  transversal  motions  which  are  governed  by  the  following  partial 
difTerential  equations  for  x  |  c  [0,f  ]: 


m 


d*tu,(f.a,)  «ti;,(f.s,)  9*ti;,(f.s,) 

+  c, - gj - EA  - - 


fiit.xi)  (2.18) 
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for  the  axial  response;  and,  for  i  =  2.3, 


dhuiit.zi)  dwi(t.zt)  d*wdt.Zi)  _  -  V  . 

m  -  +  c* - - -  +  EU  — - fiit.zi).  (2.19) 


for  the  flexural  response. 

Appropriate  boundary  conditions  for  a  "dumb  bell"  beam  structure  are 
given  by: 


EA 

8uj,(f  .0) 

8z.  * 

a2-u;,(f.0) 
dt‘  “ 

(2.20a) 

- 

EA 

8z. 

dt^  " 

(2.20b) 

for  the  axial  response,  and  for  the  flexural  response,  with  i=2,3. 

a^tiiid.o) 

dzf 

■  - 

,  aV(^o) 

dt^dzi 

0. 

(2.21a) 

r,.  flV(i.O) 

'dzf 

+ 

_  B’hv.it.O)  _ 

TTl  j  ———————  * 

0  . 

(2.21b) 

r,. 

dzf 

+ 

et^dzi 

0  . 

(2.21c) 

■  dzf 

-  - 

_  d^it.i)  _ 
*  8f2 

0. 

(2.21d) 

In  the  above  expressions  m.  mi,  mg  represent  the  distributed  mass  and  the 
two  end  masses,  respectively;  i/g.i  the  moments  of  inertia  of  the  two 
end  masses;  Cj  the  damping  coefficients;  and  EA,  EJi  the  beam  axial  and 
flexural  stiffnesses. 

The  above  set  of  partial  differential  equations,  (2.18)  and  (2.19),  may  be 
solved  using  separation  of  variables.  Accordingly,  let 


'"^(^*l)  =  SV'k(0’7i*(*i)  (2.22) 

k>i 

where  f)tt(Si)  deflnes  the  modal  shape  of  the  structure  and  )  is  defined 
as  the  modal  generalized  co-ordinate  corresponding  to  the  Ath  mode.  To 
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deflne  the  a  countably  infinite  set  of  uncoupled  ordinary  difTerential 

equations  results  through  the  projection  of  the  partial  differential  equations 
(2. 18)  and  (2. 19)  onto  the  modal  basis  : 

for  *=1,2 .  (2.23) 

where  represents  the  modal  damping  ratio,  the  frequencies  of 
undamped  free  vibration.  andpit(f )  the  modal  forces  (e.g..  see  [Clo.l]). 

For  the  models  (2.23),  linear  feedback  compensators  can  be  designed 
using  semi-infinite  optimization.  Equations  (2.21a,  b,  c,  d)  can  be  easily 
transformed  into  a  set  of  first  order  linear  differential  equations  for  use  in 
the  design  of  compensators  for  the  linear  closed  loop  control  system. 
Accordingly,  we  may  write 

*(f  .*»)  .*i)  +  G(Pb  .Pta)u(f  .*j)  (2,24a) 

y(f  .*1)  =  H(pb  ,pj|)*(f  .*,)  +  K(Pb  .Pb)u(f  .*i)  (2.24b) 

where  u  =  (ui,U8)  (for  details  of  some  previous  applications  in  beams  see 
[Tay.l]). 

(iv)  Formulation  of  Design  Constraints:  l^uige  Displacements. 

We  assume  that  the  LSS  is  to  be  designed  taking  *  specified,  large 
amplitude  changes  in  attitude  (manoeuvers)  into  account.  Although  one 
may  consider  both  fixed  and  free  time  manoeuvers,  we  simplify  exposition 
by  restricting  ourselves  to  Axed  time  manoeuvers  only.  These  manoeuvers 

can  be  specified  by  an  initial  time  f't.  a  final  time  f^t,  t  =  1,2 . *,  and 

two  linear  equations  in  the  state  vector  x  at  these  times,  of  the  form: 

L*s(f*i,s,,u,.Ue.pk  .Pb)  =  l^i  fori  =  1.2 *;  (2.25a) 

l/x(f/t,S|.u,.Ue,pb>l>b)  =  t/i  fori  =  1.2 . *.  (2.25b) 

Our  datign  variables  for  the  large  manoeuvers  are  *  pairs  of  control 
functions  Ui,«,Ue.t,  i  =  1.2,...,*.  and  the  structural  and  actuator 
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parameters  and  p^. 

During  the  large  manoeuvers  a  number  of  physical  constraints  must  be 
taken  into  account. 

(a)  Control  and  Parameter  Constraints. 

Tlje  simplest  constraints  involve  bounds  on  the  structural  and  actuator 
design  parameters  and  on  the  controls,  which  have  the  form: 

«  Pk  ^  ft  (2.26a) 

«  ft  ^  ft  (2.26b) 

l|uj.t(-.*i)IL^^^(2'i)-  j  ~  1.2:  zie/o(ft):  i  =  i.2....k  .  (2.26c) 

where  the  dependence  of  location  of  the  control  forces  is  indicated  by  the 
use  of  the  set  valued  function  /,  (the  set  /o(ft)  contains  a  finite  number  of 
points. 

(b)  Stress  Constraints. 

Next,  certain  structural  limits  must  not  be  exceeded.  For  example, 
the  maximum  values  of  the  stress  resultants  f  and  stress  couples  m  must 
be  kept  below  values  which  would  cause  yielding  or  other  forms  of  damage. 
These  result  in  inequalities  of  the  form 

-  m»(f  ,z,,u,.uz,ft,ft)<0, 

V»iC[0.1];  1  =  1,2 *,  (2.27a) 

and 

fs(<.*i)  '  f»(<  .*i.u,,ue.ft.ft)<  0, 

V*ie[0,l];  i  =  1,2 ib  ,  (2.27b) 

where  and  represent  "yield"  or  "failure"  values  and  mt  and  Are 

vectors  of  stress  couples  and  stress  resultants,  respectively,  during  the 
i—ih  manoeuver.  Both  of  these  depend  on  the  state  vector  x  and  hence  they 
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depend  on  both  the  control  inputs  and  the  structural  and  actuator  parame¬ 
ters. 

In  addition,  it  may  be  necessary  to  limit  relative  values  of  the  motion 
to  prevent  damage  to  on-board  equipment,  sensors,  or  measuring  instru¬ 
ments.  Relative  displacements  may  be  defined  as  the  difference  between 
yi(f  ,Si)  and  the  motion  of  some  "reference"  location  Jt(f.Si.ra/)  (e.g.  ^\.nj 
may  be  located  at  the  center  of  mass  of  the  undeformed  beam).  Accord¬ 
ingly,  we  define 

Ayi(f  .si.Ui.Ua.p^.Pto)  4  yi(f.2i.u,.U2.i^  .1%)  -  (2.2Ba) 

and  the  relative  translational  motions  by 

AWi(f  .Si,Ui.U2,|\  ,pb)  =  Ayi(f.zi.U|.Ue,pw  .k) 

*  <Ayi(f.Zi.U,.Uz.pk.Pi,).Ti(f.e,„/)>Ti(f.s,^/)  (2.2Bb) 

and  impose  the  constraint  that 

+  <AW*(f.S,.U,.U8.pk.Ph).AWi(f.Z,.Ui.UB.pw.Pb)><  0. 

VziC[0.l];  i  =  1,2,...,*  .  (2.2Bc) 

where  34^  is  a  limiting  motion. 

(c)  Constraints  on  Afibrations. 

Finally,  the  controls  Uj,  Ug  must  be  such  as  to  ensure  that  the  motions 
"quiet  down"  as  a  manoeuver  nears  completion.  Ihe  structural  and  actua¬ 
tor  parameters  must  facilitate  this  task.  Hence  we  get  constraints  of  the 
form 

ft(x(t,*,  u,.ug.p^.pk) )  «  0.  Vs,e[0,l],  V  f  i  s  1.2 . k.  (2.29) 

(t)  Ftarmulation  of  Design  Constraints:  Sknall  Displeeonionts. 

For  small  motion  control,  we  propose  to  build  a  tiro-degrees  of  freedom 
control  system,  as  in  Fig.  6,  with  finite  dimensional  eempensators  whose 
free  elements  form  an  additional  design  vector  .  ‘Hiare  seem  to  be  three 


id 


■saiiMnfilieaMiaH 
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main  considerations  in  small  motion  behavior:  stability,  disturbance 
suppression  and  accuracy.  In  addition,  speed  of  response  may  be  of  impor¬ 
tance.  In  the  case  of  a  finite  dimensional  system,  (see  [Pul.  ll],  attached), 
the  requirement  of  stability  can  be  expressed  either  as  a  semi-infinite  ine- 
quedity,  by  means  of  the  Extended  Nyquist  Stability  Criterion  [Pol. 6]  or  by 
imposing  constraints  on  the  eigenvalues  of  A.  Disturbance  rejection  is 
ensured  by  imposing  semi-infinite  constraints  on  the  disturbance  to  output 
transfer  function  of  the  closed  loop  system,  Hyi{ju)  over  a  critical  range  of 
frequencies,  and  so  forth. 

In  the  case  of  a  control  system  with  a  distributed  plant,  such  as 
(2.24a),  (2.24b),  the  computational  aspects  of  the  above  techniques  are  still 
to  be  worked  out  and  will  form  part  of  the  proposed  research. 

(w)  Cmnplete  Integrated  1^  Optimal  Design  Problem. 

Various  cost  functions  can  be  considered  in  the  design  of  an  LSS,  with 
the  simplest  one  corresponding  to  minimization  of  the  maximum  energy 
used  in  a  set  of  k  large  manoeuvers.  A  substantially  more  sophisticated  cost 
function  can  be  used  to  maximize  stability  robustness  with  respect  to 
modeling  errors,  etc.  In  order  to  appreciate  the  full  complexity  of  the 
optimal  design  problem,  it  is  instructive  to  state  the  model  problem 
corresponding  to  the  energy  minimization  cost: 

minimize 


I  /  ^  (2.30a) 

subject  to 


L*z(f*t,s,,ui,ue,pb  ,|]b)  =  b^t  fort  =  1.2 . k 


(2.30b) 
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l/s(t^i,Si.Ui.U2.pb.Pb)  =  tori  =  1.2 . *:  (2.30c) 

A  <  Pb  <  k  :  (2.30d) 

A.  <  Pb  «  FW  :  (2.30e) 


|uj.t(  .*,)L^O,(2,).  V8,e/.(Pb)c[0.1].  ;  =  1.2:  and  i  =  1.2 . ^2.30f) 

Wyit.Zi)  -  m«{«.2,.u,.U8.Pk,.Pb)<0. 


V2,e[0./J;  V  i  =  1.2 . *;  (2.30g) 

-  f«(t.*i.«i.U2.pb.i:b)<0. 

V2,e[0.1];  V  »  =  1.2 Jk  :  (2.30h) 

+  <Awt(<.2,.u,.Ue.pk.Pk).Aw»(1.2,.ui.ii2.pk.pb)>«  0. 

V*ie[0.i].  'r/ i  =  1.2 . it;  (2.30i) 

<'4(«(l,i|U,.Ua.pb.pb)  )  <  0.  V*ie/«(Pb).  V  t i  =  1.2 . *(2.30j) 

P^(Pfc.Pb.Pb.;«)  ^  0.  V  »g[o',m”].  j  =  1.2....m.  (2.30k) 


where  /a(Pb)  is  s  Unite  set  of  points  at  which  the  actuators  are  located,  and 
the  last  set  of  constraints  is  of  a  form  which  results  from  linear  control  sys¬ 
tem  design  (see  [Pol.  1 1].  attached). 

(vi)  Simulation  Considerations. 

The  solution  of  the  set  of  partial  difTerential  equations  introduced  in 
(iii)  and  (iv)  above,  presents  considerable  difTiculties.  even  in  the  absence  of 
the  control  or  design  aspects.  Accordingly,  it 

la  normally  neceaaary  to  utilize  numerical  methods  to  perform  the  simu¬ 
lations.  Ihe  finite  element  method  m^  be  oonveniently  used  to  develop 
an  approzimate  solution  to  the  dUIerenUal  equations  where  the  accuracy 
of  the  approzUnation  is  controlled  by  the  number  of  elements  used  in  the 
analysis,  (e.g..  see  [Zie.l]).  The  efbet  of  using  the  finite  element  method 
is  to  transform  the  problem  from  an  infinite  dimensional  one  (i.e.. 


z(  z  I ,  t ))  to  a  finite  dimensional  problem  (i.e..  z(  (z  i)( ,  t ),  where  i  denotes 
a  nodal  position  in  the  finite  element  model). 


